

REMARKS 



The present amendment makes editorial changes and corrects typographical errors in 
the specification, which includes the Abstract, in order to conform the specification to the 
requirements of United States Patent Practice. No new matter is added thereby. 

Attached hereto is a marked-up version of the changes made to the specification and 
claims by the current amendment. The attached page is captioned " Versions with Markings 
to Show Changes Made." 

In addition, the present amendment cancels original claims 1-5 in favor of new claims 
6-10. Claims 6-10 have been presented solely because the revisions by crossing out 
underlining which would have been necessary in claims 1-5 in order to present those claims 
in accordance with preferred United States Patent Practice would have been too extensive, 
and thus would have been too burdensome. The present amendment is intended for 
clarification purposes only and not for substantial reasons related to patentability pursuant to 
35 U.S.C. §§101, 102, 103 or 112. Indeed, the cancellation of claims 1-5 does not constitute 
an intent on the part of the Applicants to surrender any of the subject matter of claims 1-5. 

Early consideration on the merits is respectfully requested. 




William E. Vaughan \y 
Bell, Boyd & Lloyd LLC 
P.O. Box 1135 



(Reg. No. 39.056) 



Chicago, Illinois 60690-1135 
(312) 807-4292 
Attorneys for Applicants 
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VERSIONS WITH MARKINGS TO SHOW CHANGES MADE 



SPECIFICATION 
TITLE OF THE INVENTION 
METHOD FOR DETERMINING THE AMPLITUDE AND PHASE ANGLE OF A 
MEASURING SIGNAL CORRESPONDING TO A CURRENT OR VOLTAGE OF AN 
ELECTRICAL POWER SUPPLY NETWORK 
BACKGROUND OF THE INVENTION 

D e scription 

M e thod for D e termining th e Amplitud e and Phase Angel of a Measuring Signal 
Corr e sponding to a Current or Voltage of an Electrical Power Supply Network 

The present invention relates to a method for determining the amplitude and phase 
angle of a measuring signal corresponding to a current or a voltage on an electrical power 
supply network by using sampled values of the measuring signal, a model of the measuring 
signal containing at least a sinusoidal component being used to calculate the model amplitude 
and phase parameters of the measuring signal with the sampled values by applying a recursive 
least-squares estimation method. 

A method of this type is described in an article by M.S. Sachdev and M. Nagpal "A 
recursive least squares error algorithm for power system relaying and measurement 
applications", IEEE Trans, on Power Delivery, Vol. 6, No. 3, July 1991. In the case of this 
known method, sampled values are formed from a measuring signal corresponding to a 
current or voltage of an electrical power supply network and, by m e ans of via a linear least- 
squares estimation method, by using a sinusoidal signal model that models the measuring 
signal, are used to form a complex vector^ which specifies the amplitude and phase angle of 
the measuring signal. In this case, in a first step the real and imaginary part of the vector are,, 
in each case A determined on their own. Then, in a second step, by m e ans of via a coordinate 
transformation, the polar coordinate representation of the complex vector, that is to say the 
magnitude and phase of the vector, can be determined from the real and imaginary part of the 
vector. 

In the case of the known method, it is assumed that the frequency of the measuring 
signal is known. If this is not so, or if the frequency changes, then a separate method is 
needed in order to determine the frequency of the measuring signal. Methods are known 
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there, for example, which measure the distance of the zero transitions of the measuring signal 
and, on the basis of measuring this period, determine the frequency of the measuring signal 7 ; 
see, for example, E. Schriifer (Editor): "Lexikon MeB- und Automatisierungstechnik" 
[Measurement and Automation Encyclopedia], VDI-Verlag, 1992, p. 204. Also known is a 
5 method for frequency measurement in which the measured signal to be examined is s in each 
case^ filtered in parallel by a high-pass filter and an all-pass filter (German patent 
DE 42 1 1 946). The frequency of the measuring signal can be determined via the ratio 
between the amplitudes of the output signals from these two filters. 

The present invention is bas e d on the obj e ct of directed toward specifying a method 
10 with which all the significant variables of the measuring signal can be determined 
simultaneously and quickly. 

M SUMMARY OF THE INVENTION 

»=* 

□ In the case of a method of the type indicated at the beginning, this object is achieved, 

i^T according to the present invention, in that use is made of a model of the measuring signal 

*¥5 containing the sinusoidal component in accordance with the relationship y = A ■ sin(2 Tift + <p), 

Ly 

\ji y designating an instantaneous value of the model of the measuring signal, A the amplitude, f 
q the frequency, cp the phase angle and t the timer-b y. By using this model of the measuring 
^ signal and by using the sampled values, by moans of via a recursive nonlinear least-squares 
I y estimation method, the model frequency parameter of the measuring signal is also determined 
llO by the estimation together with the model amplitude parameter and the model phase angle 
parameter. 

Although it is known, from the book by H.J. Hermann, "Digitale Schutztechnik" 
[Digital Protection Engineering], 1997, p.p. 1 10-1 1 1, to use a recursive, nonlinear least- 
squares estimation method in protection engineering, the book does not contain any reference 

25 to the fact that by using such an estimation method, sampled values of a measuring signal can 
be used to determine the amplitude, phase angle and frequency of the measuring signal in a 
single measured- value processing process. 

A significant advantage of the method according to the present invention is, however, 
precisely that the sampled values of the measuring signal are used to determine the frequency, 

30 as well as the amplitude and phase angle, in a measured-value processing processr-and 

therefor e . Therefore, the amplitude, phase angle and frequency of the measuring signal are 
associated with the same point in time. 
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Furthermore, the patent DE 42 05 300 CI discloses a method with which the phase 
angle and the amplitude of a periodic signal can be determined via a phase-locked control 
loop (PLL (Phase-Locked Loop)). 

The use of the model for the measuring signal y = A • sin(2/rft + (p) leads to good 
5 results if the measuring signal has a purely sinusoidal waveform. If there is a DC component 
present in the measuring signal, a model of the measuring signal in accordance with the 
relationship y = A • sin(2^ft + <p) + d is advantageously used, the summand d modeling the 
DC component of the measuring signal. 

If the measuring signal is a signal whose frequency changes over time, in the case of 
1 0 such measuring signals without a DC component, a model of the measuring signal in 
accordance with the relationship 

( » \ 



I - : 
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y = A ■ sin 



^ /=0 J 



M- can advantageously can be used, and in the case of measuring signals with a DC component, a 
I a i model of the measuring signal in accordance with the relationship 



r 

y = A- sin 



3 ean advantageously can be used, f* !) designating the ith time derivative of the frequency and 

modeling a change in the frequency over time, and various orders of the time derivative of the 
frequency being taken into account by selecting the variable n. In these models, in an 

20 extension of the model specified above, y = A • sin(2;rft + (p), the frequency f is replaced by 

Z(/ ( V) 

the expression ,s=0 ; if only a 0th order derivative is taken into account, the integral 

-expression_becom_es.f (0) or f. These expanded models make it possible to determine frequency 
changes f° } over time, in addition to the variables amplitude A, phase angle (p and frequency f. 
In a further advantageous r e fin e m e nt embodiment of the method according to the 
25 present invention, the values of the amplitude A, the phase angle 9 and the frequency f 

determined by the estimation method are output as results only when the estimation error is 
less than a smallest permitted estimation error. This has the advantage that, in particular, the 
values estimated at the start of the method and afflicted by large estimation errors are not 
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output, and therefore the large estimation errors cannot have any negative consequences for a 
user of the method. 

Additional features and advantages of the present invention are described in. and will 
be apparent from, the following Detailed Description of the Invention and the Figures. 
5 BRIEF DESCRIPTION OF THE FIGURES 

In ord e r to explain the invention furth e r, Eigure 1 shows, by using a block diagram, 
the progress of an exemplary embodiment of the method according to the present invention 

Figure 2 shows waveforms of a measuring signal, of an estimated signal y s determined 
for the measuring signal by using a model y = A • sin(2;rft + <p) and of the resulting error F 7 . 
10 Figure 3 shows results of the method according to the present invention in the case of 

a purely sinusoidal measuring signal when use is made of an appropriate model of the 
measuring signal, separated in accordance with the time-based waveforms of amplitude, 
phase angle and frequency^ 

Figure 4 shows results of the method according to the present invention in the case of 
15 a measuring signal subject to an offset, that is to say containing a DC component, separated in 
accordance with the time-based waveforms of amplitude, phase angle and frequency, when 
use is made of a model of the measuring signal that does not model any DC component^and. 
Figure 5 shows results of the method according to the present invention in the case of 
ry a measuring signal subject to an offset, that is to say containing a DC component, separated in 
1^0 accordance with the time-based waveforms of amplitude, phase angle, frequency and DC 
component, when use is made of a model of the measuring signal that models a DC 
component. 

DETAILED DESCRIPTION OF THE INVENTION 
According to Figure 1, after being sampled in a sample-and-hold circuit 1 and after 
25 analog-digital conversion in an analog-digital converter 2, sampled values y m of a measuring 
signal u m are present on an input 3 of a data processing system 4. The data processing system 
4 contains a unit 5 for carrying out a recursive nonlinear least-squares estimation method; 

At the beginning of the recursive estimation method, starting values SW for the 
variables amplitude A, frequency f and phase angle cp are input into the unit 5 and are present 
30 on the output 6 of the unit 5 as estimation output values . The estimation output values 

0^ are led to a block 7, which contains a model of the measuring signal. 
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In the block 7, the estimation output values ©^ are used, taking into account the 

model of the measuring signal, to determine a starting signal value y s0 , which is output at the 
output of the block 7 and led to an input 8 of the unit 5. 

The unit 5 uses a (first) sampled value y m of the measuring signal u m and the starting 
5 signal value y s o to determine an estimation error F est j m in accordance with the relationship (3) 
recited below. This estimation error will lie above a permitted smallest estimation error. 
Therefore, in the unit 5 and starting from the estimation output values , in accordance 

with the rules for recursive nonlinear least-square estimations, new estimated values © are 

— i' ) 

formed in accordance with the relationship (1) below, are output at the output 6 of the unit 5 
H 0 and led to block 7. 



In block 7, the new estimated values © are used to determine an estimated signal 

— \') 



value y s i, which is output at the output of block 7 and is led to the input 8 of the unit 5. The 
unit 5 uses the (first two) sampled values y m and the starting signal value y s o and the 
estimated signal value y s i again to form an estimation error F es tim in accordance with the 
H5 relationship (3) recited below. This estimation error will generally also lie above the smallest 
IjJ permitted estimation error. Therefore, in the unit 5 and starting from the estimated values 
□ © on the preceding estimation, new estimated values © are again determined in 

|S| 1 (1) V ) 

accordance with the relationship (1). These steps are repeated until the estimation error F es tim 
lies below a smallest permitted estimation error. 

20 In order to determine whether the estimation error F es tim lies below a smallest 

permitted estimation error, use is made of a block 9. Via an input 10, block 9 obtains from 
the unit 5 the sampled values y m of the measuring signal u m , the starting signal value y s0 and 
the estimated signal values y s i to y S k, which are buffered in the unit 5. When the estimation 
error-Festim lies belowa-smallest pennitted estimation error, then a switching signal, which 

25 switches a switch 12 on, is output at an output 1 1 of the block 9. The switch 12 then 

forwards the estimated values © present at the output 6 of the unit 5 to an output 13 of the 

data processing system 4. The estimated values © are then output as the resulting values 

— (^) 

from the estimation method and constitute the values for amplitude, frequency and phase 
angle of the measuring signal, estimated with sufficient accuracy. After that, the estimation 
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method continues to run as described above and permits continuous determination of the 

estimated values © 

— (*) 

The unit 5 determines the estimated values © in accordance with the following 

— ) 

relationship (1). 

In relationship (1), © represents a vector which contains the estimated values of the 

— (£) 

variables to be determined, that is to say, here, the amplitude A, frequency f and phase angle 

<p, after k estimation steps; © designates a vector which results in the case of the 

— i) 

estimation after k-1 estimation steps. The matrix P (k) is a symmetrical precision matrix, as it 
is known, whose determination is presented further below using the relationship (2). The 

function ^^<*>*~ CM) ^ contains the model of the measuring signal. The vector 

h[<p J 

contains the partial derivatives of the model (k) of the measuring signal with 

respect to the parameters of the model of the measuring signal; that is to say, the derivatives 
with respect to amplitude, phase angle and frequency. The variable y (k) is the kth sampled 
value of the measuring signal. 

In order to determine the precision matrix P( k ), use is made of a relationship (2). 



E<k -" X + y T P ~ ' 7 £l *~' S 



The factor X determines the exponential weighting of preceding sampled values of the 
measuring signal At the Beginning of the method,-a starting-value_is_us.ed for the precision 
matrix P (k ). 

The estimation error F es tim after the kth estimation is determined by m e ans of via the 
relationship 
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Here, N is the number to be taken into account of the sampled values y m of the 
measuring signal u m and of the estimated signal values y s , y m i is the ith sampled value of the 
measuring signal u m and y si is the ith estimated signal value. 

The estimation error F es tim is determined by N sampled values y mi of the measuring 
signal u m and N estimated signal values y S i being evaluated. An appropriate magnitude for N 
is the quotient of the (sampling frequency of the sample- and-hold circuit 1) / (estimated 
frequency), but N ean also can be selected to be larger. During the first estimation pass, the 
value used for N is the quotient of (sampling frequency of the sample-and-hold circuit 1) / 
(starting value SW of the frequency). When, at the beginning of the estimation method, N 
sampled values y m i of the measuring signal u m , or N estimated signal values y si are not yet 
present, then only the present values are used to determine the estimation error F es tim. 

Further details relating to carrying out recursive nonlinear least-squares estimation 
methods are described in the documents J. Wede; D. Werner: "EchtzeitprozeGmodelle auf der 
Basis von Parameterschatzverfahren" [Real-time process models based on parameter 
estimation methods], automation engineering series, Vol. 214, VEB Verlag Technik, Berlin, 
1985, p.p. 30-34, p.p. 44-50, 56, 57, and Hans- Joachim Herrmann: "Digitale Schutztechnik: 
Grundlagen, Software, Ausfuhrungsbeispiele" [Digital Protection Engineering: Principles, 
Software, Design Examples], VDE Verlag 1997, p.p. 104-1 13, and A. Jurisch: "Digitale 
Impedanzmefiverfahren auf der Basis von Identifikationsmethoden" [Digital Impedance 
Measuring Methods Based on Identification Methods], Dissertation, Technische Hochschule 
Zittau, 1990, Annex 4.6.3. 

Figure 2 shows the time-based waveform of a sinusoidal measuring signal u m , whose 
amplitude, frequency and phase angle are to be determined by using the estimation method 
according to Figure 1. Also shown over time, as the curve y s , are : the estimated-signal -values- 
y s i ... ysk which are output at the output of block 7. Also shown is the waveform of the values 
of the error F which results from the difference between the measuring signal u m and the 
estimated signal values y s i ... y S k, over time t. The waveforms shown result from an 
estimation method which is carried out by using the models = A • sin(2 7rft + <p) for the 
measuring signal. It can be seen that the waveform of the curve y s is the same as the 
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waveform of the measuring signal u m after about 20 ms; the error F then assumes very small 
values. The estimation method^ therefore, permits the correct determination of the amplitude, 
the frequency and the phase angle of the measuring angle after about one period of the 
measuring signal u m . 

5 In Eigure 3, taking into account the waveforms shown in Figure 2, the estimated 

values ... 0 (jt) linked with the estimated signal values y s i ... y S k are plotted against the 

time t, separated in accordance with the variables A, cd and 2nf and cp that determine them. In 
an upper graph, Figure 3 shows the waveform of the amplitude A determined by the 
estimation method according to Figure 1, in a central graph the waveform of the variable 
10 co = 2nf that is proportional to the frequency f, and in a lower graph the waveform of the 
phase angle cp, in each case plotted against time. According to Figure 2, the result that the 
time-based waveforms have reached a steady state at their final value after about 20 ms with 
the smallest permitted estimation error, and indicate the result of the estimation method, 
y Figure 4 shows, in a manner of representation which is the same as Figure 3, the time- 

: fi 5 based waveforms of amplitude A, frequency f and phase angle cp determined in a further 
;L estimation by using the model of the measuring signal^ = A * s\n(27tft + <p) if use is made of 
sampled values of a measuring signal which contains a DC component. The waveforms of 
frequency and phase angle also contain large time-based fluctuations for times greater than 
r! 20 ms. These fluctuations point to an erroneous estimation and are a reference to the fact that 
20 the selected model is not beneficial for a measuring signal containing a DC component. 

Figure 5 shows the time-based waveforms of amplitude A, frequency f, phase angle cp 
and DC component d determined in an estimation using the model y = A • sin(2/rft + cp) + d 
for the measuring signal, if use is made of sampled values of a measuring signal which 
contains a DC component. By using the model just mentioned of the measuring signal, 
25 fluctuations in the time-based waveforms of frequency and phase angle are avoided; 

consequently, an estimation which is better with regard to the estimation eiror is carried out... 

Although the present invention has been described with reference to specific 
embodiment, those of skill in the art will recognize that changes may be made thereto without 
departing from the spirit and scope of the invention as set forth in the hereafter appended 
30 claims. 
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ABSTRACT OF THE DISCLOSURE 

Abstract 

Method for d e termining the amplitude and phas e angle of a measuring signal corresponding 
to a curr e nt or voltage of an electrical pow e r supply network 

Th e inv e ntion relates to a A method for determining the amplitude and phase angle of 
a measuring signal corresponding to a current or a voltage on an electrical power supply 
network by using sampled values of the measuring signal, a model of the measuring signal 
containing at least a sinusoidal component being used to calculate the amplitude and phase of 
the measuring signal with the sampled values by applying a recursive least-squares estimation 
method. In order to be able to determine the frequency of the measuring signal together with 
the amplitude and the phase angle, use is made of a model of the measuring signal in 
accordance with the relationship y = A ■ sin(2;rft + cp), and by using this model and the 
sampled values (ym), the determination of the frequency of the measuring signal (u m ) as well 
is carried out by m e ans of via a recursive nonlinear least-squares estimation method. By 
expanding the signal model, measuring signals with a DC component and with frequencies 
that change over time can also can be investigated. 
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